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Abstract 

It is shown that certain aspects of gravitation may be described using a relativistic 
action-at-a-distance formulation. The equations of motion of the model presented are 
invariant under Lorentz transformations and agree with the equations of Einstein's the- 
ory of General Relativity, at the first Post-Newtonian approximation, for any number 
O" 1 of interacting point masses. 

feb; PACS No. 03.30. +p, 04.25. -g, 04.25. Nx, 04.90. +e. 

Introduction 

After the discovery of the action-at-a-distance formulation of electrodynamics |T] - [TO] , 
several relativistic non- instantaneous action-at-a-distance theories have been investigated 

■ [H] - [21] • Instantaneous action-at-a-distance formulations have been studied using a vari- 
ety of approaches [22] - [SI]- For gravity, several relativistic action-at-a-distance models have 
been proposed [32] - [IS], and compared with observations [471 HE]- A major difficulty with 
most of the models proposed is their disagreement with Einstein's theory of General Relativ- 
ity (GR) in the so-called "slow motion approximation" (first Post-Newtonian approximation 

■ (1PN)), even for the simpler case of two point masses (N = 2) . 
The objective of this paper is to present a relativistic action-at-a-distance description of 

gravitational interactions for a system consisting of an arbitrary number N of point masses. 

The model presented in this paper is in agreement with GR at 1PN for an arbitrary 
number N of interacting point masses. 

Our description also agrees with GR in the so-called "fast motion approximation" for N 
point masses (at the first Post-Minkowskian order (1PM)) and it is in agreement with GR 
for the one body case (N = 1) at all orders (assuming the central mass is not spinning) if 
the Schwarzschild metric is expressed in the isotropic gauge. 

An action functional for gravity in the relativistic action-at-a-distance formulation 
In order to describe a relativistic system of N point masses interacting gravitationally, 
we consider the following action functional: 

d\d + EE ~ / / d\id\j5 (pij) Fij 

+ E E E G2m ^ mk JJJ dXidXjdX k 8 ( Pij ) 5 (p jk ) F ijk + ... (1) 
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In (JT]), rrii (i = 1, 2, N) is the mass of particle i, Aj is a Poincare invariant parameter 
labelling the events along the world line (Aj) of particle z, c is the speed of light and G the 
universal gravitational constant. 

The functions Fy F%j Ty; Tii) ? Cj); Fijk Fijk £ife> Cjfc; Tij> Tiij 7*fc> Tfc? ; 7jfc; Tfcj! C*> Cj; Cfc)> 
are invariant under Poincare transformations since they themselves are assumed to be func- 
tions of the Poincare invariants 7^, and Q. The Poincare invariants p^, £y, 7^, and 
are defined as follows [19] : 

Pij = (*i - «i) 2 (2) 
tij = {iiZj) (3) 
7« = ~ Zi)) (4) 

Q = % (5) 



We denote if = -5k The metric tensor: r\^ v = diag(+l, —1, —1, —1). 

The action functional ([ID is invariant under Lorentz transformations and does not involve 
any fields to mediate the interactions between the masses. The particles interact with each 
other directly and we assume that the interactions propagate at the speed of light c in 
vacuum. The Dirac delta functions in ([1]) account for the interactions propagating at the 
speed of light forward and backward in time. 

The action ([!]) can be written in a compact form as follows: 



S = -^rriiC J dXiQ 

i 

+ EEE- E gfc "^ 3 "^ j j r dK-dx ik l[' 

k=2 ti 12^11 ik^ii,—,ik-l '=1 



Notice that not only the two-body interactions (k = 2), but all possible k-body interac- 
tions (k = 2, N) contribute to the action. 

Without loss of generality we can assume that Fji = F^-, F k ji = F^ k and so on (^i fc i 2 ...i fc _ 1 i 1 1 



p. . . 



From (|T|), we can see that we can write the action of an individual particle i as follows: 
I' I 2G f 

I (IX; C T E "' i J d ^j S (Pij) F ij 

~xE E m i m k J J dXjdXt {5{p ij )8{pj k )F i j k + S(p jk )8(p ki )F jki + 5(p ki )5(p ij )F kij ) + ... 



Si = —rriiC 
G 2 



(7) 



The equations of motion of the relativistic particles can be derived from the action ([T]) 
(or from (171) ) using the variational principle. We find: 
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+7T7 X! £ m i m fc / / dXjdXk \S(p jk ) (-^— {8{p ij )F ijh ) - [S{pi^ 01 " 



2c 4 ^ fe ^^. J J V V V 

9 d / dF ■ ■ \ \ 

= o. 

Integrating by parts and taking into account that: 

dAi (^) dAj 7ji dAj 

we can write the equations of motion (JSJ) in the form: 

G 



■•/' 

z 



+ -J Y ">J j dX 3 5 (Pii) ( A ij + B ij*™ + 

G 2 r r 

+ — H Y m i m k J j dXjdXk5(fiij)6(j) jk ) (A? jk + B^ k i iv + Cff k z jv + Dg k z k u 

G 2 r r 

+ t^aYY m i m k J J dXjdXkSipjijSipik) (A% k + m k hv + C% k z jv + D^Zfa, 



0, 



where 



dF 



d 2 F 



13 zl 



dz ifl dz^dzif, 



Iji \ 



dF, 



ifj, t 



+- 



Iji 



~Zj Fij + 



- , ; ~ y7) ^ -r <^ ^ -r ^ 



dz 



3 

d 2 F 



in 



d F^ ^ 
dzj dzi^ j 



dz ifl dz 



'OF, 



4) 



dz 



Iji 



Fij I + 



7i< 



a 2 F, 



Iji 



dz. 



in j 



A n dF ijk d 2 F ijk Q (, a fi\ rp I dF ijk \ 

+ — l—z? F -u + (z? - z^) 3Fijk z v + £. 9Fijk + ry. . d 2F »jfc in 
+ Jji { 3 vjk + [ l j) dz] Z > + dz* + 7 « dz]dz itl z > , 



l/J, 



Jjk | / M [l\ ®Fjj k .rj 9 Fjj k , r \ f-\A\ 

\ {Zi Z3) ^4 Zk + ll3 d4dz- Zk ) ' [ } 

dz^ (15) 



jiw _ ( z i ~ z j) ( 9Fjjk _ i z i ~ rf) p \ ^ J%j ( d 2 F ijk ^ (zj — Zj) dFjj k 
%3k Jji V dz i" Tji / Tji \dzipdzjv jji dz i{ 



in 



r u = Jjkitf - Z V) ( ( z i - Z J) p lij dF ijk \ _jjk( {Z? - zj) dFy jjj d 2 F ijk 
vk llj V la 7i» 9z in ) Jkj \ Jji 9z ku iji dz ifl dz ku / 



(17) 



i J_ | i fja „v\ ®Fjik ^ t dF^ d Fji k A 

jji { 3 " h [ 1 3 ] dz) Zj + K ' 3 0%, + li3 dz]dz ilx j J 
i // „t*\ ®Fjik . v c dFjjk d 2 Fj ik \ 

+ ^ { k jik + {i ~ k) ^4 Zk + 6fc ai~ + llk ^4dz7, k ) ' ( } 

dz^ (19) 

Qixv = ( z i ~ z j) ( dFjik _ ( z i - Zj) p \ + Jiif d 2 F jik _ ( Z j - Zj) dF jik \ ^ 
Jlk Jji V diju Jji Jlk J Jji \dz ifl dz ju jji dz ipu ) ' 



Multiplying (jHJ) by ij M (and performing the summation over fx), we find that the solutions 
of the equations of motion (jHJ) must satisfy the following iV conditions (i = 1, 2, N): 
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^ ^ + ^ g ^ y [f* - ^ ) 

+ ~a12 12 m i m fc / / dXjd\ k (s(pij)6{p jk ) \F ijk - Zi-r^r 



c 



+$(Pjk)5(pki) \F jki - Zi-gTjr J + 5(pki)5(pij) \ F kij - Zi-Q^f JJ+ -J 
= 0. (22) 

Let us assume that i^...^ (A; = 2, iV) are homogeneous functions of degree two in z^, 
z ik , i.e., we assume that they satisfy the following conditions: 



■H — n-ik _ _ 2,(1 _ ii-tk _ o p /oq\ 

11 dz? --- Zi * dz? ~ Lt ^ [Z6) 

The conditions f )22|) combined with f )23|) guarantee that, for the solutions of the equations 
of motion, the expressions: 



ZG f 

Ci - —J2 m i J d ^3 S (Pij) F i3 

~ ~T 12 12 m i m k I I d\jdX k {5(pij)5(p jk )F ijk + 5{pj k )5{p k i)F jki + 5{p k i)8{p ij )F kij ) + ... 
= Ci, (24) 
are constants (which by simple scaling can be made equal to 1): 

2G 
G 2 



Q - -^J2 m i J dX Apij) F H 

rrijuik J J dX j dX k (5(p ij )5(p jk )F ijk + 5(p jk )5(p ki )F jki + 5(p k j)5{pij)F Hj ) + 



c 4 

j^i k ^ij 



(25) 



From (|23|) it immediately follows that: 



l// J /-'„ i 1 O-F; 



Using ([52]) we immediately see that the action for particle i ([7]) can be rewritten in a 
compact form as: 

Si = -m t c J dXig$%% (2T) 

where, 
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g(i) 



G x - 

^E m J 



d\jd(pi 



d 2 F 



'j 



dz?dz? 



q2 

~^iE E m i m k 

+$(Pki)8{Pij 



dXjd\ k 5(pij)8(p jk ] 



d 2 F 1 



ijk 



dz^dzf 



+ S(p jk )5(pki / 



d 2 F 



jki 



diidzf 



d 2 F, 



i- * j 



From (128 



dz?dz? 



261) it follows that: 



dzf 







dz?d% 



o. 



0. 



(28) 

(29) 
(30) 
(31) 



Using (J2 
form: 



3~T|) the equations of motion (JSJ) can also be written in a more compact 



9$% + 








(32) 



From ([TO 



From 



(33) 



, we find that the conditions ( 122]) can be simply expressed as: 

it follows that for the solutions of the equations of motion: 

rfA? = sjgdtf dtf (34) 

but also on zA 



zj and Zj (j ^ i) 



Notice that oM * s n °t a field. It depends not only on Zi, 

The main task in our formulation is to determine the functions i^-, F^^, etc, in (jH]), and 
to verify that the predictions of the theory are in agreement with observations. 

Test particles and the formulation of the action- at- a- distance model as a field theory 

Let us assume that in the limit m, — > the tensor gfW does not depend on z%. Only in 
this limit, in which rrii is a test particle, we may have a field interpretation for the metric 
tensor oW. 

Let us consider a system of iV + 1 point particles, one of them being a test particle of 
mass m and the other N particles having masses m, (i = 1, N). Let z(X) be the worldline 
of the test particle. From (l27j) we see that we can write the action for the test particle as 
follows: 



S 



-mc 



dXg 



z»z v . 



(35) 
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In fl35l) . the metric tensor g^ u depends on z (and on %i and Zi {i = 1, ...,N)), but does 
not depend on z. It can be given a field interpretation, if one desires to do so|54j. 
From (f3"2"j) we see that for a test particle the equations of motion are: 

Assuming that the matrix g^ v is invertible, these are, of course, the well known equations 
for geodesies. 

It may be possible to impose conditions on the functions Fij, Fijk, etc, in fl6]) if, for 
example, one demands that the metric tensor g^ u (which is associated with a test particle) 
obeys Einstein's field equations. Of course, there is no guarantee that this can be done at 
all orders in the Post-Minkowskian expansion, either due to the mathematical complexity 
of the equations or to the possibility that Einstein's theory of General Relativity may not 
exactly admit a dual (action-at-a-distance) formulation, or at least not one described by an 
action of the form ([6]). At the first Post-Minkowskian order it is known that = 2& — Cfj 



A possible expression for F^ 

Let us assume that the functions F^ can be expressed as follows: 



F^ = a (ey , eji) f ?, + p (ey , e i4 ) Cl , (37) 



where, 

Gm 



and 09]: 



We assume that the functions a and j3 are symmetric: 



dj = yJ&Q, (39) 
% = (40) 



a Ua , e 



" ( ( j>-<u) ■ ( 41 ) 



P(eij,e ji ) = P(e ji ,e ij ). (42) 

T/ie one-body problem 

Let us consider the case of a test particle interacting with a particle of mass M. This 
is the case N — 1 (the one-body problem). The motion of the mass M is not affected by 
the presence of the test particle. The mass M moves with constant velocity in any inertial 
reference frame. 

Let us, for simplicity, consider the inertial frame in which the mass M is at rest and 
positioned at the origin of the coordinate system. In this frame of reference the world line 
of the test particle is described by the four- vector z^ = (ct,r). From (J2"gj) and fl3Tj) we find 
the components of the metric tensor g^ u in this reference frame to be as follows: 

2GM(a(0,e) + /3(0,e)) 
(?oo = 1 ^ , (43) 

9oi = 0, (44) 
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9a = -<% ( 1 nr-^ ) • ( 45 ) 



err 



In (P-05J): 

GM 



2c 2 r 

If we choose the functions a and (3 as follows: 

[l + e) 4 - 



(46) 



4 _ (1-e) 2 
(1+6)2 



a(0,e) = i7^^> (47) 

W0= 1 ~4 e +C)4 » ( 48 ) 

one can easily check that the metric (1431 - 145]) . with a and (3 given by (147f48j) . coincides with 
the well known Schwarzschild metric of GR in isotropic form|55j. 

Since a, f3 and are Poincare invariants, we can write the functional relations: 

a(0, e j{ ) = , (49) 

At the second Post-Minkowskian order (up to terms proportional to G 2 in the metric), 
we can write: 

a(0, e^j) «a + aifji, (51) 

0(0, e if ) « fo + Piep, (52) 

where a 0) A)> «i and /?i are constants. 

The values of these constants can easily be determined by expanding ( 149]) and ( 150]) . We 
find: 

«o = 2, (53) 

A) = -1, (54) 

1 . . 

ax = --, (55) 

Pi = ~\ (56) 
The second Post-Minkowskian approximation 

Let us now consider the gravitational iV-body problem described by the action ((6]). As- 
sume that Fij are given by (137] l4T|42l . At the second Post-Minkowskian (2PM) order the 
functions a and j3 will be given by the expressions: 

~ «o + «i (e^ + e^) , (57) 
/3(ey , e^) « #, + A (ey + e*) . (58) 
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Let us consider the case where the functions can be written as: 



F%j k d\Sij ; tjii ^kii ^iki ^jk: ^kj^Cij^jkCki 

^kii ^iki ^jki ^kj)CiCjCk 
~l - c(€jjf, 6jfj, 6fcj, 6jfc, Cj'fc ; ^kj)CkiCj • (59) 

We assume that the functions a, b and c are symmetric in the indexes (ik): 

®(^kj j €jki ^iki ^-kii ^(Qj) ^jii ^kii ^iki ^jki ^kj) i (60) 

b(^kj ) ^j'fe) Cjfc, 6fcj, Gj'i, Ey) ^(^ijj ^ji; ^fci> ^ifcj ^j'fej ^kj) i (61) 

c(Cfcj') ^jfc? Cjfc; ^fci? ^j'jj ^ij) '> ^fci) ^ifej Ej'fcj ^kj) ■ (62) 

At 2PM we have: 

a ( e «i' e i*' e ^*' e ife' e jfc' e ^i) ~ a 0) (63) 

b(tij, €ji, €ki, €ik, €jk, Ejfej) ~ &(b (64) 
c(Qj') e ifc; e fcj) ~ c Cb (65) 

where a , 6 and c are constants. 

Therefore, at the second Post-Minkowskian order we can write the action, for a system 
of iV particles interacting gravitationally, as follows: 



S = - g niiC J dXiQ + g ^2 m ' mj J J dXidXjd (pij) Fij 

i i j^i C 

+ ^T 2 J J J d\dXjd\ k 5 {pij) 5 {p jk ) F ijk , (66) 



where, 



i j^i k^i,j 



^ = (ao + «i (ey + eji)) ej + (/?o + A {tij + Cj (67) 



^iife — aoiij^jk^ki + bodCjCk + c o^fciCi (68) 

At the second Post-Minkowskian approximation there is no need to consider the functions 
Fi x ___i k for & > 3 since the terms associated with these functions in the action ([6]) give 
contributions only at the (A; — l)-Post-Minkowskian order. 

At the second Post-Minkowskian order (2PM) the equations of motion are: 



+ — Y,Y, m i m k I d\jd\ k 5{pij)8(p jk )A\£ 



+ i^jY.J2 m i m kJ j ' d\jd\ k 5{pji)8{p ik )Af^ 

= 0. (69) 
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Substituting flST l 168]) into ffTTI - IT3]) and (1X4 ] [T8|) we find (in this approximation): 

i0)/< V"' '■')(,. .<■'! - J > 2 / & J | % ) 



.4 



(i)m 



+ 



+ 



2 ( 



z? - Z M 



CiVji 

2z[_ 

l 

CjVji 



^4 (| + 1) ft - (i - 1) M + 3 ftc. 



, Vij Vji , 



B 



ij 



a 



(0)/ii/ 



CjVji 



(«o4 + AC?) + 



Z; — Zj 



CjVji 



CjVji 



(aodjZj + PoCjZ- 



,VijC, 



VjiCj 



+2^\ (a (Trey + ^r) + 2 a^) , 



^ijk 







1 



~ ( + — C/c ) ) (aoiijCjkiki + froCiCjCfc + c o£fc;Cj 



7j< 

0— + 6j - — f + — Ck ) ) OCfc^o 



Iji \ Iji 



Ikj 



Ikj 



Z; 



Jij lij 



Ijk, 



H — ( [ Cj— ~ ( tjk + — Ck I I Cjk£ki a O — CiCjCkbo — CkiCj c O 

Iji \\ Iji Ikj V Ikj 



lij 



lij 



H — — I Cj— + £ij — I + — Cfe ) I { a oijkdj + 2c £fciC/) 

7ii V 7ii 7fci V Ikj ) J 
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7. 



("ui/,C/.'s/v + b (j(i(k + coCfcjO 



7i 4 



+ 2if (0^ + <^ + £ + £) (OCA. + 
+i i ( [CrY + C*nr + — J 6y&*ao - — (^oOCiCfc + coCfcjC 

V \ 7ji 7fci 7fci / 7ji 

\\ Iji 1h IjiJ In v 



(75) 



From ( IBTl l68l) and (|28l) it follows that, at the second Post-Minkowskian approximation, 
the metric tensor associated to the particle i (with non-negligible mass mi) is given by the 
formula: 



5^ — ^ ^ 



^ mj y d\ j 5(p ij ) [a z j ^z ju + PoCjTfpv] 



G 2 r S(p--) 



G 2 m, 



J2 m i / dX 



\Vji\ 



+ 



Z%ii \Z%v Zjv) 



hj\ 



Zip, \Zifi Zjp) 



lij 







lij 



Zip, \Zip Zjfi ) 



lij 



G 2 f f 
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+5(p jk )5(p ki ) ( -aoCife ( z knZjv + ZjpZkv) + b Q(kV^ + caCkZjpZj, 



+5(pki)5(pij) (y a o€jk {zkfiZjy + z j)jL z ku ) + (boCkCj + coCffc) 



Z%u \Ziu Zj v 



lij 



(76) 



For the case of a test particle (the mass of which can be neglected) in the presence of N 
particles with non- negligible masses mj (i = 1, N) the above expression simplifies to the 
following: 



2G 

9pv Vfiv 



J2 m i d\i5((z - Zij 2 ) laoz^Ziv + /3 z 2 rj^ 
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c 4 ^ ^ !(*(*-*)) 



dXidXj 



8((z - Zi) 2 )5{{zi - z,-) 2 ) f ~a (ijij) + z jfl z iv ) + M, ? Zjty» + coz 2 z jfl z 

+5{{zi - Zjf)5((z - Zj) 2 ) (~a (ziZj) (z ifjL z ju + z jfl z iu ) + b^z 2 ^ + c z 2 z^z^ 
+5((z - ^-) 2 )(5((2; - Zi) 2 ) (~a (ziZj) {z ifl z jiy + z jtl z iv ) + (b z 2 z 2 + c (ij%) 2 ) 



(77) 

JTie /irsi Post-Newtonian approximation 

The equations of motion flo^]) involve multiple times. The force acting on mass i depends 
on the state of motion of particle i at time t and, to account for the time needed for the 
transmission of the interactions, on the states of motion of the remaining N — 1 particles at 



the past and future times tl (j ^ i, s = — , +) and also on t ( ^' s> (k ^ i,j, s = — , +). 

Using Taylor series expansions involving the particles' present motions at time t, one can 
rewrite the equations fl69|) using just the one time variable £[56j. We use series expansions 

up to terms of second order (\) (first Post-Newtonian approximation (1PN)). 
From the definition ([5]) it follows that: 



d\j 



rdtil-f" 



Cdt, ( 1 - ^ 



(7f 



(79) 



The Dirac delta function can be expressed as follows [TO] : 

/ 



5(c 2 (t 



~ (n - rj) 2 



In f lHUj) . tj (s = — , +) are the two roots of the equation: 

c\t-t,) 2 -(f i (t)-f J (t J )) 2 = 



1 

2c 



S(tj - tf~ h 



+ 



Ah+h 

"3 °j 



5{tj - tf 



V 



R ret _ V » 1 J R adv + V g J 



+) 



(80) 



and, 



(*-«n 

(f +, -0 
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(81) 

(82) 
(83) 



Rrf = f. _ ^.") (84) 

fljf = rl - rf ,+) (85) 

£ — n 1 ' •* is the time it takes for a signal to travel forward in time at the speed of light 
from particle j to particle i. 

tj^ — t is the time it takes for a signal to travel backward in time at the speed of light 
from particle j to particle i. 

To terms of second order we can write: 



n-^' « m + + Vj^^ - a^, (86) 



n-rf^^r^-v/^+v^-a^ (87) 

^-f, (88) 



From (I81ti89p . we find (to terms of second order): 



ffret 

ij 



i ^ (^) 2 (njfljO N ^ 

2c 2c / 



1 ij 



1 - -i 



i 

(i,+)2\ 2 



Jj>a<fo _|_ V 13 3 J 



2c 2 2c 2 



(91) 



From the definitions HI EJ) and ( 1391 1301) . it is not difficult to see that, to terms of 
second order, we can write: 



Cy + 2c 2 + 2c 2 c 2 ' 19 j 



(* -) w r . . ( ! + (^-) 2 , (fkM 



,r-« , (93) 



\2 



, r „_ r + ^ + + , (95) 
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Vij ~ r v{ l+ c c + 2c 2 + 2c 2 2c 2 J' W 

In (1561 - 1961) . = r*i — r} is the relative position of particle i with respect to particle 

j, Hij = — , t4 is the velocity of particle i, and Vj, Sj the velocity and the acceleration of 

particle j. All these quantities are given at time t. 

From (I25p and ([57)1 . at the first Post-Minkowskian order , for the solutions of the equations 
of motion we obtain: 

r „ c/ N> / £ 2 



C» = 1 + E m j J dXjSip^Cj ( a -f + A) J ( 97 ) 



Now, substituting (179|80p and (I90|l91|l92p into ( 197]) . to terms of second order, for the 
solutions of the equations of motion we can write: 

Ci«l + ^(ao + A))E— • (98) 

Substituting (1751- 196]) and (}9~5p into (1691- 175]) . we find the equations of motion to terms 
of second order (in \) (first Post-Newtonian approximation): 



at +G {a + Po) 2^ —riij + —di H — H — Oj 2^ (3a - Po) E — a i 

~ TTi ( Q o + A)) £ —nijijUjaj) + % E -T™iJ ( a o ( v i + v i ~ 2 (^)) ~ | («o + Po) (n^) 2 

^ . ^ jji . 

+ — (Po - ao) % E ~t ((jHjfii) ~ (jkjVj)) + — E ~2^i {2a {n ij v i ) - (a - (3 ) (n^Vj)) 



+ — ^ ((a + Po) 2 + «i + Pi) E + — ( 2 («o + Po) 2 + «i + ft) £ -r 
c i# '''/ c jV< 

£r 2 ^-^ Tfljm^ / 1 / o \2 \ 1// « \2 \ l 

+ tLL 2 n »i — I 2 ( a o + Po) + «o + o + c J H ^(ao + Po) + «o + o + c J 

= 0. (99) 

Complete agreement with the equations of motion of General Relativity [571 158] , at the 
first Post-Newtonian order, is achieved if: 



a = 2, 


(100) 


A) = -i, 


(101) 


«! + /?!= -2, 


(102) 


ao + &o + c = -2. 


(103) 
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Conclusions 

We have obtained Lorentz invariant equations of motion describing the gravitational 
interactions of a system consisting of iV point masses. The equations are derived explicitly 
from a Lorentz invariant action. Contrary to General Relativity, which is a field theory, the 
model presented here is a relativistic action-at-a-distance description (the interactions are 
not mediated by a field). We have shown that the equations of motion for N point masses 
agree with those of General Relativity at the first Post-Newtonian approximation (1PN). 
Agreement with General Relativity for the N body problem at orders beyond 1.5PN has not 
been established. The model presented is in agreement with General Relativity for the one- 
body case, at all orders. At the first Post-Minkowskian approximation our model reduces to 
the model of Havas and Golberg [501 EE], which is known to be in agreement with General 
Relativity in this approximation. Due to this agreement, gravitational radiation effects in 
our model begin to appear at the 2.5PN order (^-) [5H1 ED]- 
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